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Abstract 

The generators and commutation relations are calculated explicitly for 
higher symmetry algebras of a class of hyperbolic Euler-Lagrange systems 
of Liouville type (in particular, for 2D Toda chains associated with semi- 
simple complex Lie algebras). 
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Introduction. We give a description of the generators and relations in higher 
symmetry algebras for a class of Darboux-integrable hyperbolic Euler-Lagrange 
systems of Liouville type [U [51 [3] . There exist many non-equivalent definitions 
of this type of PDEs [H El Hj ; we investigate the systems f l that admit as many 
first integrals of the characteristic equations Dy{w) = and Dx{w) = on i?l as 
there are unknown functions. The 2D Toda chains u^y = exp(i^tt) associated 
with semi-simple complex Lie algebras are the most well studied example of 
such equations [U [31 13 [HI [7] • The systems of this class are known to possess 
higher symmetries = 0(0) that depend on free functional parameters (f) = 
*(0i(a;, [w]), . . ., (jirix, [w])) and belong to the image of matrix total differential 
operators □ (linear operators in total derivatives) [H [3 [9l [10] . The existence of 
such operators □ for Liouville-type systems was observed in [H [3] and [31 [TO] , 
where the importance of the linearizations i'^^ of the first integrals w in the 
construction of □ was revealed. In the paper [H] we proved that the additional 
assumption for be Euler-Lagrange strengthens known results and even makes 
the description of □ explicit, see formula ([3]) below. 

In this paper we establish the transformation rules for the operators □ under 
unrelated reparametrizations of the coordinates in their domains and images. 
We show that, under natural assumptions on the geometry of iSl, the images 
of these operators are closed with respect to the commutation, whence the 
Lie algebra structure on their domains appears. We calculate the brackets on 
the domains explicitly, which yields, by the push forward of the Lie algebra 
structure, the commutation relations in the symmetry algebras sym^L- To do 
this, wc introduce auxiliary Hamiltonian operators which have the same domain 
as □. 
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Remark. We do not assume the presence of a symmetry x <-> y in £l- We 
work with 'the x-half of the algebra sym^L related to the first integrals G 
kerDj,|g ; the reasonings hold for the respective 'y-half of sym^L, and the two 
subalgebras commute between each other. For the Euler-Lagrange systems £l 
at hand, the integrals S keiDxl^ are not used in the proofs, unlike in [10] 
for arbitrary Liouville-type systems. 

The full list of assumptions on the systems £l and their integrals is given 
in our main Theorem [HI see also Remark [5] on p. [HI However, the reasonings 
in section 1 hold under less restrictive conditions. In particular, the number 
of first integrals w^, . . ., w"^ for the characteristic equation on can be less 
than the number of the unknowns u^, . . ., u™ in £i^. In that case, the auxiliary 
(r X r)-matrix operators Ak defined in Q become smaller in size but remain 
Hamiltonian (see [S] for the second Poisson structure for KdV provided by the 
2D Toda chains with a unique integral). 

The paper is organized as follows. First we define the operators □ that de- 
termine symmetry generators for the systems £l and introduce auxiliary Hamil- 
tonian operators. Here we re-derive the higher Poisson structures for the Drin- 
fel'd-Sokolov hierarchies TT] on 2D Toda chains related to semi-simple complex 
Lie algebras; an example is given for the A2-Toda chain. Then in section 2 we 
establish the commutation closure for images of the operators □ and calculate 
the structural relations in the algebras sym^L! an illustration is given for the 
Kaup-Boussinesq equation. Finally, in section 3 we discuss some properties of 
the operators that yield symmetries of non-Lagrangian Liouville-type systems. 

All notions and constructions from geometry of PDE are standard fT2j [131 
[14] . We follow the notation of O [151 IE] ■ This paper develops further the 
concept of [S]- 

1. Symmetry generators for £l- 

Definition. A Liouville-type system f is a system {u^y — f{u,Uj:,Uy;x,y)} 
of m hyperbolic equations upon u = (u^, . . . , u™) which admits nontrivial first 
integrals 

, . . . , £ kei Dy\^; ui^, . . . , w'^ € ker_D^|^, < r,f < m, 

for the linear first order characteristic equations Dy | ^ (w* ) = and \ ^ (w^) = 
that hold by virtue (=) of £. 

Example 1. In [2j it was proved that the 2D Toda chains fS] u^.^ = exp{K^jU^) 
related to semi-simple complex Lie algebras with the Cartan matrices K admit 
maximal (r = f = m) sets of the integrals. Various methods for reconstruction 
of w;* , for these exponential-nonlinear Toda chains were proposed in [H [H [7] . 
The differential orders (after a shift by —1) of the integrals w^, . . ., w"^ w.r.t. 
u are equal to the exponents of the corresponding semi-simple complex Lie 
algebras of rank r, which follows from [2", p. 21]. 

For instance, in the sequel we consider the Euler-Lagrange 2D Toda system 
^Toda associated with the simple Lie algebra s[3(C), see [UlSll^, 

^^Toda = = exp(2u - w), = exp(-lt + 2w), K=[_\^D]. (1) 
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The integrals of respective orders 2 and 3 for system ^ are (e.g., see [T7]) 

= Uxx + Vxx -ul + UxVx- vl and w"^ = u^xx - '2uxUxx + UxVxx +ulvx- Uxvl . 

The generators 1^9 = □(</)(a;, [w])) of higher symmetry algebras for Liouville- 
type equations are given by matrix total differential operators □, see [31 [3]. For 
Euler-Lagrange Liouville-type systems £1, — {F = E(£) — 0}, see [H1[H1[TH], the 
existence of certain factorizations for at least a part of symmetries is rigorous 
and can be readily seen as follows. For integrals w such that Dy(w) — V(F) 
vanishes on the differential ideal {F = 0}°° by virtue of an operator V, and 

for any I{x, [w]), the generating section Tpi = V* o o (1) for a 

conservation law J I dx solves the equations ^e(£)(V'/) = on £l, see [HI [131 
[T3]. The Helmholtz condition (e{c) — ^e(£) ^'^^ ^^^'^ linearization (the Frechet 
derivative) implies that the vector 

(^[0] = [V* o (4"))*] {d^ix, H)) e ker€E(£)|£^ (2) 

is a symmetry of iSl for any — (£j"''')*(l) — 'Ew{I dx). A standard homological 
reasoning (see [131 Ch. 5] or §7.8]) shows that formula ^ yields symmetries 
of the system £1^ even if sections do not belong to the image of the variational 
derivative E.^, w.r.t. w. 

In this section we recall the construction of operators □ that determine 
symmetries for a class of Euler-Lagrange Liouville-type systems. We suppose 
that the integrals w are minimal, meaning / € kerZ3y|g implies / = I{x, [w]). 

Proposition 1 ([9]). Let k be an invertible symmetric constant real (to x to)- 
matrix. Suppose that C = J L dxdy with the density L — 
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Hi,{u;x,y) is the Lagrangian of a Liouville-type equation £l = {E(£) = 0}. 
Let m = dL/duy be the momenta and w{m) — {w^, . . . ,w^') be the minimal 
set of integrals for £1^ that belong to the kernel of F>y\^ . Then the adjoint 
linearization 

□ = (4"T (3) 

of the integrals w.r.t. the momenta yields Noether symmetries ipc of £1^: 

ipc = n{5'Hl5w) for any H = j H{x, [w]) dx. (4) 
Corollary 2. Under the assumptions and notation of Proposition[Tl the section 

(^ = n(0(x,H)) (5) 

is a symmetry of the Liouville-type equation £l for any r-tuplc = *(0i, . . . ,(j)r). 

Proof. Consider the jet bundle J°° (^) over the fibre bundle ^ : M'' x M ^ K with 
the base R 3 a; and the fibres W with coordinates w^, . . . , w''. By Proposi- 
tion[Tl the statement is valid for any in the image of the variational derivative 
E^. Obviously, its image contains all variational covectors whose compo- 
nents 0i(x) e C°°(M) are functions on the base of the new bundle ^. The 
prolongation of the substitution w = w[m[u]\ : J°°{n) — » r(^) converts the 
components of sections to smooth differential functions in u (which denotes 
the set of fibre coordinates in the bundle tt over the same base R 3 x). Now 
recall that □ is an operator in total derivatives Dx whose action on differential 
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functions f[u] is defined by joo(s)(£'x(/)) := ^(joo(s)(/)) through the restric- 
tions joo(s)(/) of / onto the jets joo(s) of sections u = s{x). Hence we obtain 
(j)i{x) = 4>i[x, [w[xn[s{x)]]\) , and the assertion foUows. □ 

Remark. This proof combines taking the variational derivatives with respect to 
w on one jet space with calculating the total derivatives of differential functions 
in u on the other jet space over the same base. Whenever the two bundles 
coincide, the reasoning amounts to the definition of D^- Then it is called the 
substitution principle (|13|, see a detailed discussion in [14j). 

Theorem 3. Under differential reparametrizations w = lilw] and u — u[u\ of 
the coordinates , . . . and u^, . . . , u™ in the infinite jet bundles over ^ and 
n that specify its domain and image, respectively, the operator □ is transformed 
according to the formula 

a ^ a = o D o (e^i"K* , (6) 

It I w=w\u\ \ I 

U — U\VL\ 

Proof. The transformation ip = of the velocities is obvious. Under dif- 

ferential reparametrizations w = w[w] of the integrals, the sections <p — STC/Svu 
are transformed by = (^Ib"')*!^); thence □ becomes well defined on imE^,. 
Namely, it maps variational covectors for the fibre bundle ^ to evolutionary 
derivations in the jet space over the other fibre bundle tt. Repeating the rea- 
soning used in the proof of Corollary[51 we establish the transformation rule ^ 
for □ on the entire domain. □ 

In Theorem [3] we showed that sections in the domain of the operator □ are 
transformed by the same rule as the arguments of Hamiltonian operators. There 
is indeed a deep reason for that. 

The integrals ^[m] of Euler-Lagrange Liouville-type systems £l determine 
the Miura substitutions from commutative modified KdV-type Hamiltonian hi- 
erarchies 25 of Noether symmetries for to completely integrable KdV-type 
hierarchies 21 of higher symmetries of the multi-component wave equations 
^0 = {sxy = 0}, see below. A natural example is given by the potential modified 
KdV equation ut = —\uxxx + u^, which is transformed to the KdV equation 
Wt = —\wxxx + "iwWx hy w = — Uxx- This example was analysed in de- 
tail in [S]. The method for generating relevant differential substitutions by the 
integrals w of Liouville-type systems was discovered in [T^, see [S] for discus- 
sion. This fact was used in [5Dj for a classification of the first-order differential 
substitutions. 

The hierarchies 21 and 58 share the Hamiltonians 7ii[m] = Tii [^[m]] through 
the Miura substitution ^[m]. The Hamiltonian structures for the Magri schemes 
of 21 and 05 are correlated by the operators □, which map cosymmetries (pi for the 
hierarchy 21 to symmetries </?i+i of the modified hierarchy *B. We stress that the 
first Hamiltonian operator Bi = (^m'')* for 55 originates from the differential 
constraint m = dL/duy upon the coordinates u and the momenta m for £i^. 
Using the explicit form of the 'junior' operator Bi and the differential-functional 
closure in \w\ for the velocities of the integrals, see [3], we realize the classical 
scheme for generating higher Poisson structures via Miura's substitutions |21| . 
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Lemma 4. Introduce the linear differential operator 

ife = □* o Si o □ (7) 

that maps variational covectors for the jet bundle J°°(^) over ^ to evolutionary 

vector fields on it, At: V{p ®c-(R) C^{J°°{C))-r^ r(e) ®c-(R) C^{J°°{C))- 
The operator ([7]) is Hamiltonian and determineslj a Poisson structure for the 
KdV-type hierarchy 21. The coefficients of Ak are differential functions in w. 

Proof. By construction, the Poisson bracket {Til, 7^2}^^ = (Eu,Hi, Afc(E^H2)) 
satisfies the equality 

{HiN,H2N}a, = {nAw[m]],H2[wM] (8) 

Therefore the left-hand side of ^ is bi-linear, skew-symmetric, and satisfies 
the Jacobi identity. Fourth, it measures the velocity of the integrals w along a 
Noether symmetry of £1^. Since evolutionary derivations are permutable with 
the total derivative Z?y, the velocity {Hi, 0.2} lies in kerZ?j,|^ and hence its 
density is a differential function of the minimal integrals w. □ 

The multi-component wave equation £0 = {s^y = 0}, whose symmetries 

contain the hierarchy 2t and such that Ai = (^1^')* encodes the differential 
constraint between the coordinates s and momenta w for £0 , is chosen such that 
the first structure Ai = A^^ for 21 factors the higher Hamiltonian structure for 
«8. Hence Bfe. = □ o Ai o □*, where k' = k'{D, (iw^)*) > 2. 

Example 2. Consider the Euler-Lagrange 2D Toda system ([T]). The density L 
of its Lagrangian is 

L = — i((2u2; — Vx) ■ Uy -t- {2vx — Ux) ■ Vy^ — exp(2M — v) — exp(2w — u). 

Therefore we introduce the momenta :— 2ux~-Vx and := 2vx-'Ux, whence 
we express the integrals as w = ^[m]. All symmetries (up to x <-> y) of ([T]) are 
of the form ip = □(c/)(a;, [w^], [w^])), where (p = '((/'i,(/'2) is a pair of arbitrary 
functions and the (2 x 2)-niatrix total differential operator is 
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Ux + Dx -jDl - UxDx - \ul - luxVx + + \ux 

Vx + Dx + luxx - \vxx - + luxVx + \vl 

The entries of the arising Hamiltonian operator A2 — \\Aij, 1 < i, j < 2|| are [ 
3 , o„„in I „„i 



All = 2Dl + 2w^Dx + wl 
A12 = ~Dl ~ w^Dl + {iw^ - 2wl) ■ Dx + {2wl - u;^), 



A21 = Dt+ w^Dl + iw'^Dx 
A22 = -iDl - \w'Dl - 2wlDl + (2wl - 2wlx - l{w' f) ■ Dx 
+ U^wlx - 2wlxx - 2w^wl). 

The shift w'^ 1-^ + X oi the second integral yields the 'junior' Hamiltonian 
operatoJl A^^^ := ^ (A2) — (3]^^ ^^"^ ), which is compatible with A2. 



^In most cases, this is one of the higher structures for 21, which is indicated by the subscript 
k = fc(n, m) > 2. The choice of the 'junior' operator Ai for 21 is discussed in what follows. 
^The analogous operator A^^"^ = ^|^_Q(^fc)i where + fi, is not Hamiltonian. 
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The pair {A\ , A2) is the bi-Hamiltonian structure for the Boussinesq equa- 
tion w] = 2wl ~ wl^, w1 = -Iwl^^ - Iw^wl + wl^. The symmetry = 
{A2 o 6 / 6w){^J Ax) starts the second sequence of Hamiltonian flows in the 
Boussinesq hierarchy 2t. The modified Boussinesq hierarchy 58 shares the two 
sequences of Hamihonians with 21 by the Miura substitution w — w[m] with 
m = Tn[M]. Namely, for any Hamiltonian 7i[w], the flows Ur = (5?i[m]/(5m, 
rrir — —5T-i\m[u\\ 1 5u belong to !8. The velocities Ur constitute the commuta- 
tive subalgebra of Noether symmetries for the 2D Toda chain ([T]) . 



2. Commutation relations in symi?L- In this section we prove the commu- 
tation closure for the images of operators ^ by using the well-known analogous 
property of the auxiliary Hamiltonian operators ([7|) . At once, we describe the 
relations in the symmetry algebra generated by ([5]) for Sl- 

First, consider a linear total differential operator A whose arguments 4>{x^ [w\) ~ 
'((/>!, . . . , (j)r) are the variational covectors for the infinite jet bundle over ^. As- 
sume that the image of A in the Lie algebra of evolutionary vector fields is 
closed w.r.t. the commutation: [imA, imA] C imA. By the Leibnitz rule, two 
sets of summands appear in the bracket of fields A{(l)'), A{(l)") that belong to 
the image of A: 

[Ai^'), Ai^)] = A{dAi^,0') - 9 AW) i^')) + (dAm (^) i^") - 9 AW) (^) i'P'))- 

In the first summand we have used the permutability of evolutionary derivations 
and total derivatives. The second summand hits the image of A by construction. 

The commutator [ , ] | j^^^ ^ induces a Lie algebra structure [ , ] a in the quotient 
ri(^^) of the domain of A by its kernel: 

[A{<t>'),A{<t>")\ = ^([0', rU), 0', r e m.). (9a) 



This bracket, which is defined up to ker^, equals 



= - dAW)i^') + (9b) 

It contains the two standard summands and the skew-symmetric bilinear bracket 
{{,}}a. 

Lemma 5 ([131 El)- The image of a Hamiltonian operator A = ||X]t ^t^{Xj i"^])' 
Dr\\ is closed w.r.t. the commutation. The k-th component (1 < fc < r) of the 
arising bracket {{ , }}^ on the domain of A is calculated by the formula 



\a\>Oi=l |r|>Oj = l 



The coefficients of the bilinear terms in the bracket {{ , }}^ are differential func- 
tions of the variables w. 

Now we pass from the Hamiltonian operators (O to the operators □ that 
have the same domain as Ak but take values in a different Lie algebra. Here is 
our main result. 
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Theorem 6. Let the following condition^ he satisfied on an open dense subset 
of the Euler -Lagrange system £l = {uxy = f{u;x,y)^ of Liouville type {see 
Proposition [T]): 

• the constant symmetric real matrix k in the kynetic term of the Lagrangian 
density C be invertible; 

• the linearization i^^^ — \\df^ /du^\\ = f'{u;x,y) of the right-hand side 
in i?L be an invertible matrix; 

• there be as many integrals w*(a;, [m]) G kerZ3j,|g^ as there are unknowns u^; 

• the integrals w be minimal: $ G kcrU^I^ implies $ — (x, [w]); 

• the integrals w be differential -functional independent^ meaning that $(a;, 
[w[m]]) = implies <i> = 0; 

• the {r X m)-matrix A = ||9t(;*/9tn;^^^^ ||, where d{i) := orda; is the differ- 
ential order of the z-th integral w^[m\, be invertible. 

Then the following statements hold: 

1. the image of the operator ^ is closed w.r.t. the commutation of symme- 
tries ip — \w\)) e sym^L; 

2. the bracket {{ , }}□ arising on the domain of the operator □ satisfies the 
equality 

{{</'','/'"}}□ = {{0',0"}h,, (11) 

3. the coefficients of the bilinear terms in the bracket {{ , }}□ are differential 
functions of the integrals w. 

Remark 1. The first assumption of the theorem implies that the system is 
determined, normal, and ^-normal (see section 3). We emphasize that is the 
only system of equations imposed upon the sections u — s(x,y) G r(7r). 

Our second statement means that the ambiguity (up to ker^fe) in the choice 
of a representative from the equivalence class {{(f>', <j)"}}^^ in the right-hand side 

of (jlip amounts to the choice of an element from ker □ C ker Ak ■ We prove the 
equality of the kernels for all ^([w]) G ^{^tt)- This implies that commutation 
relations in sym , which are determined by the Lie algebra structure (j9bp on 
^(Ctt), are obtained explicitly via (fTO|) for Ak- 

Being a corollary of Lemma [5] and the first two, our third statement is, at 
the same time, a special case of Proposition [7] (see below). 

Proof of Theorem ^ We notice first that symmetries ([5]) are independent of u 
and of Uy, Uyy, .... Hence this is also true for the commutator ip = [if' , if"] G 
sym^L of two such symmetries ip' = D(^(f>' {x,[w])) and (p" = 0(0" (a;, [w])) , 
because the Lie bracket is a local bi-differential operator. 

^The list can be not minimal such that it is easier to verify each requirement. 

*The non-existence of a nontrivial and hence of its nonzero linearization ^'j""' = f'jj"' o 

is equivalent to (V o with V = £^'"') V = 0. This nondegeneracy requirement is 
dual to H12| l. see below. 



7 



The factorization ([7]) and Lemma [5] provide the diagram 

□ 

The commutator (p = [ip',(p"] determines the velocity $(a;, [wfj of the integrals 
that equalfH d[^r ^^//^{w) = [Ak{4i'), Ak{4>")]- Since the image of the Hamilto- 
nian operator is closed under commutation, we obtain the equivalence class 
cj)(x, [w]) = [4>':4'"]^^ of sections such that <f> — Ak{4>) — 9n(0)(w)- By con- 
struction of Ak , the commutator of ip' and ip" belongs to the set of symmetries 
□ ((/)). This proves the first statement of the theorem. 

However, there may be many such ip = □(</>) that induce the same velocity 
w = <f>. Since kerD C ker A^, then, in principle, the equivalence class [0', 4'"]^^ 
may contain elements that do not belong to [0', 0"]n. 

We claim that all representatives of the equivalence class [(/>', 0"]^^, determine 
a unique symmetry p of ^l- Therefore this p = is the commutator 

[□(</)'),□(</)")] of the two symmetries, because the image of <j) under □ must 
contain it. It suffices to prove the uniqueness of the trivial solution p for the 
linear homogeneous equation £w\p) = (^^^1 °^m^){v) = 0- 

There exist three ways to obtain the zero velocity "I>(a;, [w]) = of the 
integrals ?i;[m[ii]] along p G sym^L- The first possibility is that the integrals be 
differentially dependent, which is excluded by the assumption of the theorem. 
Second, the intermediate equation nD^^p) = 0, det k ^ 0, may have nontrivial 
solutions only if some shifts p — const are symmetries of £h- However, the 
determining equation [D^Dy — /'(m; x, y)) (const) = on £l then exprimes 
the linear dependence between the differentials of its right-hand sides. This 
contradicts another initial assumption. 

Therefore the proof of the second statement is reduced to the uniqueness 
problem for the zero solution ip — Q oi the linear homogeneous equation ^^^^j (V') = 

0. Since is the only systenj^ of differential relations that are imposed upon 
the sections u = s{x, y) S r(7r), the zero in the right-hand side of £'i^\ip) = is 
achieved identically w.r.t. [m] (otherwise it would overdetermine £i^). 

There remain two situations when a velocity tp of the momenta m = —^ku^ 
makes w — ^ zero. One reason is obvious: it is the use of the equation i?l — 
{uxy = /}• Indeed, if the 'time' along the flow rh = '0 is the variable y, 

1. e., my = —SHi^/Su, then we have d^\w) = Dy(w[m]) = on ^l- But the 
presence of u in the list if arguments of / excludes such solutions from further 
consideration. 

^The ^-normality of implies that ifi is its symmetry whenever the velocity dip{w) of the 
minimal integrals lies in kerDj^l^^, see H16|l in section 3. 

^The hyperbolic system is formally integrable [12j : its infinite prolongation exists 
and there is an epimorphism £^ — > 
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Without loss of generality we assume that the integral w''[m] has the high- 
est differential order: d{r) > d{i) for all i < r. Let us calculate the velocities 
of the non-minimal integrals (w'Y :— oi^^^ '''^^\'w^). Using the permutabil- 
ity [Da;, 9^'"''] = of evolutionary derivations with total derivatives, from the 
identities dll^\w^) = ■ ■ ■ = d^\w^) = we deduce that d'^\w') = 0. It is 
readily seen that the linearization of the new integrals is of the form — 
A • Dx^^'^ + 0{d{r) — 1), where the matrix A is invertible by our initial assump- 
tion. Multiplying the new equation i^^^ij-^) = by A~^, we obtain by induction 
that [tti]) does not depend on [m], hence ■0 = V'(a^)- Consequently, the 
admissible sections ip that solve the intermediate equation tpi^) = —^fiDx{(p) 
also depend on x only: ip = p{x). However, such sections, whenever nonzero, 
can not b^ll symmetries of the hyperbolic system Uxy — f{u]x,y) = due to 
the nondegeneracy det/'(M) ^ 0. In this notation, the "symmetry" ip{x) must 
satisfy the determining equation (^D^ o Dy — f'{u)^ = on £l- The first 
summand vanishes because Dy(x) = 0. Thus we obtain f'{u) ■ p{x) — 0, where 
the linearization matrix f'{u) is invertible. Hence 1^9(2;) = 0. This completes 
the proof. □ 

Remark 2. In the proof, we arrived at the linear ODE (l^'' (^p{x)) = that 
holds simultaneously for all sections s G r(7r), although nonzero solutions ip{x) 
do not contribute to the symmetry algebra. This is possible, first, if there is a 
total differential operator V such that o V = 0. (For instance, the identity 
0) (-D )(^(^)) = holds for all h{x).) To avoid this, it is necessary 

• to require the nondegeneracy of the linearization of the integrals: 

^5ii°V = =^ V = 0. (12) 



In the adjoint form, V* o □ = V* = 0, equation (fT2| exprimes the 
absence of linear differential relationfH between components of symme- 
tries p = □(^(x, [w])). 

This property is dual to the nondegeneracy Vo^^' = V = that originates 
from the differential-functional independence ^{x, [w]) = ^ <!> = via V = 
i"^, see footnote m on p. [71 

Finally, let the section s{x,y) G Sol be a solution of the Darboux-integrable 
Liouville-type system Sj^. Taking the restriction — i^'' I . , , of the lineariza- 
tion operator onto the jet of s, we obtain the ordinary differential equation 
£^(V'(a;)) — 0. For each s, the linear space 0{s) of its solutions is finite di- 
mensional. (For example, its dimension is equal to the sum of the exponents 
of a semi-simple complex Lie algebra if is the associated 2D Toda chain.) 
Therefore 



the requirement 

seSol£LCr(7r) 



''This is an immediate, point-by generalization of the fact (see above) that ifi = const 7^ 
is not a symmetry of £l- 

^Thence the nondegeneracy H12|l is analogous to the notion of normal differential equa- 
tions in the analysis of their formal integrability, see section 3. 
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in combination with (jl2p . permits to eHminate the excessive freedom in the 
choice of solutions ^{x) of the equation £^\')p) = 0. 

Theorem [S] is iUustrated for semi-simple complex Lie algebras of rank two 
in [16j . where the Hamiltonian operators Ai and Ak are constructed for the 
corresponding Drinfel'd-Sokolov hierarchies [H] and the commutation relations 
in symfL are calculated for the 2D Toda chains u^y = exp{Ku). 

Example 3 (The modified Kaup-Boussinesq equation). Consider an Euler- 
Lagrange extension of the scalar Liouville equation [15 , 

A,y = -^Aexp{-\B), B,y^^exp{-\B). (13) 

Denote the momenta hy a — ^B^ and b = \Ax. The minimal integrals of 
'3' 



system ([T5|) are wi = — jo^ — ax and — ab -\- 2bx such that Dy{wi) = on 



T3|) . i = 1,2. Hence the operator 



4 Bx Dx 2 

^Bx-2D, 



determines (Noether, see (jl])) symmetries of ([T3l) . The bracket {{ , }}□ induced 
in the inverse image of □ is 

where i/j = * (ip^ , tp"^) and x — *(x^:X^); we use upper indices for convenience. 
Consider a symmetry of (11; 



At = \AxAxx + \ {\Al - 1) Bx, Bt = -2A,,, + \AxBI - \AxBxx- (14) 

Let us choose an equivalent pair of integrals u — W2, v = wi + jW2- The 
evolution of u and v along (fT4|) equals 

Ut^UUx+Vx, Vt = {uv)x +Uxxx- (15) 

This is the Kaup-Boussinesq system, and (|14p is the potential twice-modified 
Kaup-Boussinesq equation, see [22] . The right hand side of the integrable sys- 
tem ([H)) belongs to the image of the adjoint linearization □ = (^(u t^))*- The 
operator □ factors the third Hamiltonian structure = □* o (^[j^'^'^'')* ° Q 
for (fT5|) : we have fc = 3 and 



A 



KB 



uDx + ^Ux Dl + {\u^ +v)Dx + \iu^ + 2v)x 

^ + {\u^ + v) Dx + ^Vx ^{2uDl + 3uxDl + {Suxx + 2uv)Dx + Uxxx + {uv)x) 

By Theorem m the bracket {{ , }}^ is equal to {{ . }}^kb, which is given by 
formula pU]) . We obtain 

where Vi = -i dI+vD^) ^nd V2 = -i (^"^ J^^J. The operator Ai = 

{ D ^0°) ^^^^ Hamiltonian structure for (fT5|) : its inverse Ai = A^^ factors 

the second Hamiltonian structure B2 = O o Ai o O for (fH]) . 
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3. Non-Lagrangian Liouville-type systems. Let Z = \F = 0} be a Li- 

ouville-type system; now it may not be Euler-Lagrange. Let a column w G 
kerDy|^ be composed by minimal integrals for thence Dy{w) ~ V(F) for 
some operator V. By construction of the Liouville-type systems £, there are 
no differential relations (syzygies) between the hyperbolic equations {F^ = 0} 
in them: A(F) = implies A = 0. For the same reason, the systems £ are 
normal [TH [T3]: A o ip = on £ also requires A = 0. Consequently, an 
evolutionary vector field 9^ is a symmetry of a Liouville-type system £ if and 
only if it preserves the integrals, 

Dy{d^iw)) = a^(V)(F) + V(9^(F)) = V(£fM) on £. (16) 

Consider the operator equation 

Dyoi^^ = Vo£j. on 5. 

If, hypothetically, a total differential operator □ such that 

4"^ o □ e CDiff (ker Dy\^-> ker Dy\^ (17) 

were constructed, then it would assign symmetries Lp ~ n(</>) of the Liouville- 
type system £ to arbitrary r-tuples </>(a;, [w\) of the integrals, see ([5]). 

The recent paper [TU] contains an algorithm for construction of operator 
solutions □ for the equation in total derivatives 

4' o □ = l,„x„ . Dl mod CDiff<fe {]^crDy\^ ^ kcri?y| J. (18) 

Most remarkably, the truncation from below for the sequence of coefficients of 
lower order derivatives in □ is a consequence of the presence of a complete set 
of the integrals w £ kerDj^j^ for £. However, the minimal integrals w must 
be 'spoilt' by differentiating w.r.t. x a suitable number of times. Consequently, 
instead of the Hamiltonian operator Ak — 4"* oD, see ([7]), one obtains the r.h.s. 
of (fT5|) . Likewise, the images of operators constructed in piT do not always span 
the entire ^-components of the Lie algebras symf , and the images are generally 
not closed under the commutation. Moreover, the transformation rules in the 
domains of □ under reparametrizations w[w] of the integrals remain unspecified 
for non-Lagrangian Liouville-type systems. 

Proposition 7. If the image of a solution □ of the operator equation pT]) for 
a Liouville-type system £ is closed under the commutation, then all coefficients 
of the bracket {{ , }}□ on its domain, see ([9]), belong to ker Dy\^. 

Proof. By assumption, we have that (^Dy o o □)([(/)',(/)"]□) = for all 
0', (j>"{x, [w] ). This equals 

= (Z?,o4)on)(9a(^,)(0")-an(^„)(0')+{{0','^"}}n) = (4"^ oD) (i?, ({{</>', </,"}}□)), 

because the underlined composition satisfies ^T7\ . Clearly, Dy{(p') and Dy{ip") 
vanish on £ for arbitrary 0', (p" . For the same reason, not only the whole bracket 
{{(/)',(/>"}}□, but each particular coefhcient standing at the bilinear terms in it 
lies in kei Dy\^. □ 
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Example 4. Consider the parametric extension of the scalar Liouville equation, 



£{e) = {u^y = exp(2w) • ^Jl + Ae'^ul], e e K. (19) 

This equation is ambient w.r.t. the hierarchy of Gardner's deformation of the po- 
tential modified KdV equation, see [E]. The contraction 14 =U{e, [u{e)]) from 
(|T9l) to the non-extended equation Uxy = exp{2U) is U — u + ^ arcsinh(2£Wj,); 
it determines the third order integral for (fT9|) using the integral w = — Uxx 
at £ = 0. However, the regularized (at e = 0) integral of order two for is 



w 



(1 - Vl + 4^)/2e' + u,,/v/r+4^. (20) 



The second integral for (fT9| is w = Uyy — — ■ exp(4u) e 'kei Dx\^f^^y The 
operators □ = Uj, + and 

□ = i(l + 4e2«2 _ 2e^ux.) ■Dx+Ux + Ae'ul - 2s\,,x + ^^''^i ^^^^ 

assign symmetries ip = □((/)(x, [w])) and^ = □(0(j/, [w])) of (fT^ to its integrals. 

The images of both operators □ and □ are Lie subalgebras in sym£(e). The 
bracket {{p, q}}-^ = Pyq — pqy for □ is familiar [31 E]- The bracket induced in 
the domain of □ has the following form: for any arguments p, q, we have 

{{p, q}}a = ■ {pxxqx - Pxqxx) - 2e^ • {pxxxq - pqxxx) 

- I2e^ ■ (Se^uluxx - ^s'^ulu^xx + 'i^^Uxul^ + ^u^u^x ~ u^xx) 

X [l + 8e^u^ + IQe'^ui - 2e^Uxx - Se^uluxx] • {pxxq - pqxx) 
+ (l + 288e'^M;^ - 2'&d,e\luxx + 2Se^ul - IGe^Uxx - 288£^ u^UxxUxxx 

- 96e^ul^ + m2e^°ul° + 2\e^u\^,^ ^ 2Ae^u^x + MOSe^u^ + 3328£®m^ 

- l&Se^'^u/^xUxxUx - "iMe^Uixuiuxx - 2'iMe^u^x'U'xxUxxx + iMe^u^xxUxUxxx 

- 'imSe^^uluxxUxxx + 16e'*wL " 5632£8m^^^^ _ i<^2Qe^Uxxui + 3328e*M;^u2^ 
+ 512e^ul^ul + 384e^°'u;^u^^^ - 960e^°ul^ul - A8e\xUxxx - iG72e^°uluxxx 
+ 'iQ72e^°ulxui - 2-iMe^uluxxx - blQe^uluxxx + 288£^M42:U^ + 384£*m^'u^^ 

+ 6144£^°U^^M^ - 6144£^°U:ra;W^ + 1152£*U4a:M;^ + 1536£^°M4:.M^ + 192£^uLa;W 

+ 240£^u;j;^ + Ibi&e^^ul^uluxxx - A&e^u^xUxx) 



X [1 + 96£^u^ + 2b&£^ul -f 256£^u^ + 4£4"L - 48£^u^u^^ -|- 'i2e^ul^ul 



- Ae'^Uxx - 2b&e''uluxx + M£''ulu':^x " I92e''uxxul + 16£^<J • [pxq - pqx 

The two underlined units correspond to the bracket Pxq — pqx on the domain 
of the operator [I\ —Ux ^ ^^x that provides symmetries of the Liouville equa- 
tion Uxy = exp(2W) at £ = 0. In agreement with Lemma [7l the non-constant 
coefficients of bilinear terms Pxxq — pqxx and Pxq — pqx in {{p, <?}}□ belong to 
kei Dy\^^^y It is remarkable that, since the entire construction p^H^Tt contains 
formal power series u(s) in e, so are these two rational functions: an attempt to 
express their dependence on [w] leads to formal series with unbounded growth 
of the differential orders of its coefficients. 
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Discussion. The matrix operators □ = (□*'-', 1 < « < 1 < j < r) given 
by ^ are generalizations of tensors of type (2,0) in the geometry of infinite 
jet bundles. We define the operators by using the two unrelated groups of 
differential reparametrizations for the coordinates in the domains and images, 
respectively. Furthermore, the operators □ for the Liouville-type systems 
generalize the theory of Hamiltonian structures as follows: they map variational 
covectors for one equation (we recall that sym£0 D 21) to symmetries of the 
other system (such that sym^L ^ S)- 

Unlike in [Sllini, we do not attempt to solve equation p7|l upon □. On the 
contrary, we define the operators ([3]) by a geometric reasoning. Thence, first, we 
obtain the Hamiltonian operators Ak = £w o □ for the KdV-type hierarchies on 
Euler-Lagrange systems of Liouville type [IHIS] and, second, we prove that the 
images of the operators □ are involutive. In other words, we describe a direct 
algorithm aimed at constructing completely integrable equations. 

Formulas ^ and ^ prescribe the differential order of Ak- Estimates for 
the orders of the integrals w for the 2D Toda chains associated with semi- 
simple complex Lie algebras g are known from [5] , see Example [1] and were 
reformulated in [H . The upper bound, that the numbers ord^; — 1 are not 
greater than the exponents for g, is proved by verifying (via Schur polynomials) 
Serre's relations {a,dYi)~^ i'^^{Yj) — 0, i ^ j, for the generators 

of the characteristic Lie algebras (see [J dl [7] and also [16 ), and by using 
Frobenius theorem. The fact that the vector fields Yi coincide with the Chevalley 
generators fi of the semi-simple Lie algebra q is important here. The same 
estimate from below follows from the absence of relations other than Serre's for 
the generators Yi. This was established in [2; p. 21] for the root systems A„ and 
D„ by listing the linear independent nonzero iterated commutators. 
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